Bipolar thermoelectric effect in a serially coupled quantum dot system 
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The Seebeck coefficient (S) of a serially coupled quantum dot (SCQD) junction system is theo- 
retically studied via a two-level Anderson model. A change of sign in S with respect to temperature 
is found, which arises from the competition between tunneling currents due to electrons and holes 
(i.e, bipolar tunneling effect). The change of sign in S implies that one can vary the equilibrium 
temperature to produce thermoelectric current in either the forward or reverse direction, leading to 
a bipolar thermoelectric effect. For the case of two parallel SCQDs, we also observe the oscillatory 
behavior of S with respect to temperature. 



Owing to energy and environment issues, it has be- 
come important to consider novel applications related 
to the thermal properties of materials. Many consider- 
able studies have been devoted to seeking efficient ther- 
moelectric materials because there exist potential appli- 
cations of solid state thermal devices such as coolers 
and power generators. A quantum dot-(QD) based 
thermal device was also predicted to have more pro- 
nounced enhancement in energy conversion.^) Recently, 
some theoretical efforts have focused on the thermoelec- 
tric effects in nanostructure junctions, ^''^^■^^ however, not 
many works have paid attention to the thermoelectric 
effects of a serially coupled quantum dot (SCQD) junc- 
tion, which exhibits features of current rectification due 
to spin blockade, negative differential conductance, non- 
thermal broadening of electrical conductance, and coher- 
ent tunneling (for identical QDs) in the Coulomb block- 
ade regime.^''-' Our recent work has described these ob- 
served phenomena in a unified theory. Based on our 
previous work, we find that the Seebeck coefficient of SC- 
QDs exhibits a behavior of sign change with respect to 
temperature arising from electron Coulomb interactions. 
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where the numerators denote probability factors for var- 
ious charge configurations, and they are pi = (1 — 
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Ni^s{l - N,^„ - Nj^s + cj), pe = Ni^aiNj,a - Cj), pr = 
Ni^a{Nj^a — Cj), and — Ni^g-cj. The denominators for 
the eight configurations are (i) Hi — figfij 
dots empty, (ii) 112 = (m^ - Ut^j){nj -Uj) 
A empty and dot B filled by one electron with spin a, 
(iii) Ha = {fif — Uij){fij — Ujj) — tj^ with dot A empty 
and dot B filled by one electron with spin a, (iv) 114 = 
{fii—2Uij){fij—Uj — Uj^i)—t^ with dot A is empty and dot 
B filled by two electrons, (v) lis = {^g — Ue){fij — Ujj) — tl 
with dot B empty and dot A filled by one electron with 
spin a, (vi) Xlg = {fie-Ui~Ue,j){pj -Uj -Uj^i) -t^ with 
both dots filled by one electron with spin a, H-^ = (fi£ — 
Ui — Ue.j){fJ.j — '2Uj,i) — with dot A filled by one electron 
with spin a and dot B filled by one electron with spin cr, 



Using the Keldysh-Green's function technique, we and (viu) lis = (/i£ ~Ue,- 2Ue,j){p,j - Uj - 2Uj^i) - 



can express (up to the second order in the interdot cou- 
pling, tc) the tunneling current through a serially coupled 
QDs connected to metallic electrodes (shown in the inset 
of Fig. 1) 



J = derie)[fUe) 



(1) 



where T{e) = TLie)TR{e){Ai2 + An) 12 is the trans- 
mission factor. r£=i_fl(e) denote the tunnel rate from 
the left electrode to dot A and the right electrode to 
dot B. /L(fl)(e) = \l\ei^-i^i^(R))I^BTnR) _^ denotes the 
Fermi distribution function for the left (right) electrode. 
The chemical potential difference between these two elec- 
trodes is related to pL — fJ-R = eAV. Tl{r} denotes 
the equilibrium temperature of the left (right) electrode, 
e and h denote the electron charge and Plank's con- 
stant, respectively. For simplicity, we consider wide-band 
limit that is Ti{e) — Tg. Aij denotes the spectral den- 
sity, which can be calculated by one particle off diagonal 
Green's function. In the atomic limit, we have 



with dot A filled by one electron with spin a and dot 
B filled by two electrons, iig = e — Eg + iTg. The no- 
tations Ei, Ui, and Uij denote, respectively, the energy 
levels of dots, intradot Coulomb interactions, and inter- 
dot Coulomb interactions, tc denotes the electron hop- 
ping strength between two dots. 

The probability factor for all channels of Eq. ([2]) are 
determined by the thermally averaged onc-particlc oc- 
cupation number and two-particle correlation functions, 
which can be obtained by solving the lesser Green's 
functions. We have N^^^ = ^(l/"") / defe{e)ImG^^^{e), 
and C£ = — (I/tt) / defi{e)IniG\ ^(e), where the retarded 
Green functions Gj;^(e) and G^^(e) are, respectively, 
given by 
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Occupation numbers of Eqs. (3) and (4) should be solved 
self-consistently. Note that in the absence of tc the ex- 
pressions of Eqs. (3) and (4) can also be found in our 
previous works. ^^'^^ In the linear response regime, Eq. 
(2) can be rewritten as 



J = CiiAV + £i2AT, 



(5) 



where AT = Tl—Tr is the temperature difference across 
the junction. CoefScients in Eq. (5) are given by 
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Here T{e) and /(e) = l/[e('-®^)/'=«'^ + 1] are eval- 
uated at thermal equilibrium. If the system is in an 
open circuit, the electrochemical potential (AV) will be 
established in response to a temperature gradient; this 
electrochemical potential is known as the Seebeck volt- 
age. The Seebeck coefficient is defined as S* = AV/ AT = 
—C12/C11, where £11 denotes the electrical conductance. 

Using the following physical parameters: Ue = Uo = 
SOTo, U12 = lOFo, and Tl = Tr = ITq, where the 
average tunneling rate Fq has been used as a conve- 
nient energy unit, we numerically calculate the thermo- 
electric coefficients £11 and £12- Figure 1 shows the 
electrical conductance Ge and Seebeck coefficient (S) 
as a function of temperature at nonzero orbital offset 
(^1 - E2 = AE 0). For the case of E2 = Ep - Uq 
and AE = U2 — U12, the electrons injected by a small 
bias can only tunnel through the the spin singlet state 



of the SCQD. This is the so-called "spin blockade" ef- 
fect of the SCQD. ".15,18,19) rpj^^ electrical conductance 
is suppressed when E2 deviates from the resonance con- 
dition with AE = U2 — U12, and we observe that there 
is a zero-crossing temperature Tq for the Seebeck coeffi- 
cient, i.e, S{T) = at T = Tq. A zero Seebeck coefficient 
5 (To) indicates that the current arising from the temper- 
ature gradient can be self-consistently balanced without 
electrochemical potential. The negative S indicates that 
electron carriers of the left (hot) electrode diffuse into 
the right (cold) electrode via the resonant channels above 
Ep, the negative AV is built up to reach the condition 
of J = at open circuit [see eq. (5)]. For example, the 
curve of Ei = Ep — IOFq and E2 = Ep — 30Fo has a neg- 
ative Seebeck coefficient. Note that the resonant channel 
El + U12 — E2 + U2 [see the inset of Fig. 1(b)] has a 
zero contribution in S. On the other hand, the Seebeck 
coefficient is positive when holes of the right electrode 
diffuse into the left electrode via the resonant channels 
below Ep. Here, we define the imoccupied states below 
Ep as holes. Consequently, the change in sign of S is at- 
tributed to the competition between tunneling currents 
due to electrons and holes. To further clarify the mecha- 
nism of S in sign change with respect to temperature, we 
consider the case of zero orbital offset for simplicity. 

Figure 2(a) shows the Seebeck coefficient as a function 
of temperature for various electron Coulomb interactions 
with Ee — Eq — Ep — IOFq. In the noninteracting case 
{Ue = Uej = 0), S is always positive (see dash-dotted 
curve). This is because the tunneling process is domi- 
nated by holes of the hot electrode (left electrode) dif- 
fusing into the cold electrode (right electrode) through a 
level Eq below Ep in the weak tc limit. To reveal the 
mechanism of sign change in the Seebeck coefficient for 
the interacting case (as shown by solid curves in Fig. 
2), we analyze the contributions associated with differ- 
ent poles described in eq. (2). In the weak tc limit, there 
are six poles associated with the left dot for the spec- 
tral density: e = Eq, e = Eq + U12, e = Eq + 2Ui2, 
e = Eo + Uo, e = E0 + U0 + U12, and e £^0 + C^o + 2C/i2- 
In addition, we also find that only four channels (pi, p^, 
Pq, and Ps) have high probability weighting. They cor- 
respond, respectively, to the resonant channels e — Eq, 
€ = E0 + U12, € = E0 + U0 + U12, and € = Eo + Uo + 2Ux2- 
When U\2 = IOFq, the ps channel does not contribute to 
S because its pole location is aligned with the Fermi level, 
i.e, Eq + U12 = Ep. The strengths of pi, pg, and pg for 
U12 = lOFo as functions of temperature are shown in Fig. 
2(b). The hole contribution is given by pi, which leads to 
a positive contribution to S. The electron contribution 
is governed by pg and ps, which correspond to resonant 
channels with energy above Ep, providing a negative con- 
tribution to the Seebeck coefficient. Consequently, the 
change in sign of S results from the interplay between 
the competition of electron and hole flows. The behavior 
of the Seebeck coefficient near the zero-crossing temper- 
ature To is linear. When the equilibrium temperature is 
away from Tq, the sign change in the Seebeck coefficient 
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indicates that one can produce thermoelectric current in 
either the forward or reverse direction, leading to a bipo- 
lar thermoelectric effect. On the basis of the closed form 
solutions of transmission factors and Green's functions, 
we can solve £11 and £12 in terms of polygamma func- 
tions to find To accurately. We find that Tq is mainly 
dominated by in the atomic limit {tc/Uc ^ 1), and Tq 
increases with increasing Ui, at fixed E^, but is insensi- 
tive to the variation in Ui2- These results imply that the 
SCQD may have a stable Tq under small fluctuations of 
the QD size, and a high value of Tq can be achieved for 
small size QDs . 

To achieve a thermal device with a high density of 
charge and heat currents, we need to consider a higher 
SCQD density. Consequently, the proximity effect be- 
tween SCQDs on the Seebeck coefficient should be inves- 
tigated. For simplicity, we employ the case shown in the 
inset of Fig. 3 as an example. The detailed expression of 
transmission factor for two parallel SCQDs can be found 
in ref. 15, where we investigated the charge ratchet effect 
on current rectification. Figure 3 shows the Seebeck co- 
efficient as a function of temperature. The curves S\{T), 
S2{T), 53 (T), and 54(T) correspond respectively to the 
intcrdot Coulomb interaction U = 0,4,6, and SFq. In 
Fig. 3(a), we consider Ue = SOTq, Ei = Ep — lOFo, and 
U12 = C/34 = lOFo. In Fig. 3(b), we consider = 6OF0, 
Ei ^ Ep - 2OF0, and U12 ^ U34 = SOFq. We observe 
that To is pushed toward a lower temperature with in- 
creasing U. In Fig. 3(a), S4{T) (blue line) is negative 
in the entire temperature regime. Such a behavior in- 



dicates that the number of resonant channels involving 
electron Coulomb interactions above Ep increases with 
increasing U and they dominate electron carrier trans- 
port. In Fig. 3(b), we obtain UbTq = 27Fo in the curve 
of 5i(r). Compared with the black line in Fig. 3(a), 
To is enhanced. This is attributed to the increase in in- 
tradot Coulomb interactions U^- In Particular, there are 
two zero-crossing temperatures in the S2{T) and 5*3 (T) 
curves. They are fcsTo = 6F0 and fc^To = 22Fo in the 
S2{T) curve, and /cbTq = 5Fo and fc^To = IQFq in the 
53 (T) curve. The oscillatory behavior of the Seebeck co- 
efficient in the 5*2 (T) and ^3 (T) curves is observed. There 
is only one zero-crossing temperature, ksTo = 14.5Fo, in 
the S4{T) curve. 

In this study, we find that the sign of electrochemical 
potential can be tuned by selecting the equilibrium tem- 
perature for a given temperature gradient. This implies 
that a temperature-controlled bipolar thermoelectric de- 
vice can be achieved. For the two parallel SCQDs, the 
oscillatory behavior of the Seebeck coefficient with re- 
spect to temperature is also observed. 
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Figure Captions 

Fig. 1. The electrical conductance Gp and Scebcck co- 
efficient as a function of temperature for various values of 
E2 (from Ef - lOFo to Ep - SSFq) with Ei=Ef- lOFo, 
tc = O-lFo and T l — V n = IFq. Insets shown in Fig. 
1(a) and 1(b) illustrate, respectively, the SQCD con- 
nected to the metallic electrodes and the band diagram 
corresponding to Fig. 1(a). 

Fig. 2. Seebeck coefficient as a function of temperature 
&t El = E2 = Ep — IOFq for different electron Coulomb 



interactions. Other parameters are the same as those of 
Fig. 1. Diagram (b) shows the probability of resonant 
channels in the case of U12 = IOFq. 

Fig. 3. Seebeck coefficient as functions of temperature 
for various values of interdot Coulomb interaction C/13 = 
U2A = U and Ui4 = U24 = U/2. Diagrams (a) and (b) 
consider two different sets of physical parameters U( = 
30Fo, Ei = £;F-10ro, U12 = C/34 = lOFo and Ue = 6OF0, 
Ei = Ef- 2OF0, Ui2 = C/34 = 2OF0, respectively. 
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